Strain-Induced Coupling of Spin Current to Nanomechanical Oscillations 
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We propose a setup which allows to couple the electron spin degree of freedom to the mechani- 
cal motions of a nanomechanical system not involving any of the ferromagnetic components. The 
proposed method employs the strain-induced spin-orbit interaction of electrons in narrow gap semi- 
conductors. We have shown how this method can be used for detection and manipulation of the 
spin flow through a suspended rod in a nanomechanical device. 
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An ability to control the spin transport in semicon- 
ductors is a key problem to be solved towards implemen- 
tation of semiconductor spintronics into quantum infor- 
mation processing Q, 0, y] Many methods have been 
proposed to achieve control of the electron spin degree 
of freedom using magnetic materials, external magnetic 
fields and optical excitation (for a review see Ref. 0). 
Other promising ideas involve the intrinsic spin-orbit in- 
teraction (SOI) in narrow gap semiconductors to manip- 
ulate the spin by means of electric fields Q and electric 
gates jM 0, • Recently, some of these ideas have been 
experimentally confirmed 0,0]. 

In semiconductors the spin-orbit effect appears as an 
interaction of the electron spin with an effective magnetic 
field whose direction and magnitude depend on the elec- 
tron momentum. A specific form of this dependence is 
determined by the crystal symmetry, as well as by the 
symmetry of the potential energy profile in heterostruc- 
tures. In strained semiconductors new components of 
the effective magnetic field appear due to violation of 
the local crystal symmetry [l0| . The effect of the strain- 
induced SOI on spin transport was spectacularly demon- 
strated by Kato et. al. in their Faraday rotation exper- 
iment ||. An interesting property of the strain-induced 
SOI is that the strain can be associated with mechan- 
ical motion of the solid, in particular, with oscillations 
in nanomechanical systems (NMS), in such a way mak- 
ing possible the spin-orbit coupling of the electron spin 
to nanomechanical oscillations. At the same time a big 
progress in fabricating various NMS |nj allows one to 
reach the required parameter range to observe subtle ef- 
fects produced by such a coupling. 

In this Letter we will consider NMS in the form of a 
suspended beam with a doped semiconductor film epi- 
taxially grown on its surface (see Fig. An analysis 
of the SOI in this system shows that the flexural and 
torsion vibrational modes couple most efficiently to the 
electron spin. As a simple example, we will focus on the 
torsion mode. The strain associated with torsion pro- 
duces the spin-orbit field which is linear with respect to 
the electron momentum and is directed perpendicular to 



it. This field varies in time and space according to re- 
spective variations of the torsion strain. Due to the linear 
dependence on the momentum, the SOI looks precisely as 
interaction with the spin dependent electromagnetic vec- 
tor potential. An immediate result of this analogy is that 
the time-dependent torsion gives rise to a motive force on 
electrons. Such a force, however, acts in different direc- 
tions on particles with oppositely oriented spins, inducing 
thus the spin current in the electron gas. The physics of 
this phenomenon is very similar to the spin current gener- 
ation under time-dependent Rashba SOI, where the time 
dependence of the SOI coupling parameter is provided 
by the gate voltage variations |fj. In the present work 
we will focus, however, on the inverse effect. Due to the 
SOI coupling, the spin current flowing through the beam 
is expected to create a mechanical torsion. The torque 
effect on NMS due to spin flow has been previously pre- 
dicted by Mohanty et. al. ^3] for a different physical 




FIG. 1: Schematic illustration of electromechanical spin cur- 
rent detector, containing a suspended semiconductor-metal 
(S-M) rectangular rod atop an insulating substrate (blue). A 
spin current is injected from the left semiconductor reservoir 
(yellow) and then diffuses toward the metallic film (green). 
While passing through the semiconductor film, the spin cur- 
rent induces torque shown by the black arrow. 
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realization, where the torque has been created by spin 
flips at the nonmagnetic-ferromagnetic interface. They 
also suggested an experimental setup to measure such a 
small torque. As it will be shown below, the torque due 
to the strain induced SOI can be large enough to be mea- 
sured using the experimental setup proposed in Ref. 01 ■ 
Besides this method, other sensitive techniques for dis- 
placement measurements can be employed 

The system under consideration is a rectangular beam 
of the total length L t , width 6, and thickness c. The 
coordinate axes are chosen as shown in Fig. 1. The semi- 
conductor film with the thickness c/2 occupies the length 
L of the beam. The rest part contains a metal film. It 
can also include some additional elements for detection 
of the torque, for example in Ref. Here we will 

consider an example when the spin current is created by 
diffusion of the spin polarization from the left contact 
in Fig. H Therefore, there is no electric current flow 
through NMS. The spin polarization diffuses towards the 
metal film which, due to its relatively high conduction, 
can play an important role as a reservoir for the spin 
polarization relaxation. 

We start from the strain-induced SOI 0] described by 
the Hamiltonian 

H S oi = a[a x (u zx k z -u xy ky)+a y (u X yk x -Uy Z k z ) 

~\~ (7 2, i^yz^y ~~ ^zx^x Uyy U zz ) 

+ ayky {u zz - u xx ) + a z k z (u xx - Uyy)] , (1) 

where Uij arc elements of the strain tensor, oi stand for 
Pauli matrices and k% denote components of the electron 
wave vector. In the narrow gap semiconductors the pa- 
rameter (3 is usually much smaller than a floj |. There- 
fore, the term proportional to (3 will be omitted below. 
Besides the strain-induced -ffsoi, the total SOI Hamil- 
tonian also includes the strain independent interaction 
Hso2- Because of submicron cross-section dimensions of 
the doped semiconductor film, i?s o2 will be determined 
by the bulk Dresselhaus term |l5| . 



H S0 2 = 5Y,\z iin \h{k?j-kl) 



(2) 



This interaction, in the range of doping concentrations 
10 17 cm -3 and higher, provides the main mechanism for 
spin relaxation in bulk materials [lflj . 

Since the S-M rod with total length L t ^> b and c, 
the major contribution to the strain comes from flexural 
and torsion motions of the rod [J^ ■ Within the isotropic 
elastic model the flexural motions are represented by the 
diagonal elements u xx and u yy [l4| which do not enter 
into the first square brackets of Eq. (JTJ. On the other 
hand, due to the crystal anisotropy effects, the u xy com- 
ponents are not zero for such sort of motion and could 
contribute to Eq. Q. We will consider, however, the 
simplest example of torsion motions of the rod within an 



isotropic clastic model. In this case the strain can be 
represented as |l4j 
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where t{x) = 88 /dx stands for the rate of torsion deter- 
mined by the torsion angle 6. The function \ depends 
only on z and y and is uniquely determined by the rod 
cross-section geometry. 

The next step is to derive from the one-particle inter- 
action Eq. m a Hamiltonian which describes a coupling 
of the spin current to the strain. The electron system 
carrying the spin current can be described by a density 
matrix p. In the framework of the perturbation theory 
the leading correction to the electron energy due to the 
SOI induced strain can be obtained by averaging -ffsoi 
with p. In the semiclassical approximation such a pro- 
cedure can be represented as averaging over the classi- 
cal phase space with the Boltzmann distribution func- 
tion -Fk(r)- This function is a 2x2 matrix in the spinor 
space. One can also define the spin distribution func- 
tion f£(r) = (l/2)Tr[F k (r)cr 1 ]. It is normalized in such 
a way that the local spin polarization J"(r) = ^ k Pk(r) 
We notice that, due to electron confinement in y and z 
directions, the averages of i^soi containing k y and k z 
turn to zero. Assuming that electron distribution is uni- 
form within the cross-section of the semiconductor film 
one thus obtains, from Eqs. JI} and 10, the SOI energy 
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This expression can be further simplified ta king into ac- 
count that x turns to zero on a free surface [1J|. Hence, 
in the example under consideration x = on the top and 
side surfaces of the doped semiconductor film. Conse- 
quently, the second term in Eq. vanishes after inte- 
gration over y. Now Eq. Q can be expressed in terms of 
the spin current J y (x) which is the flux in x-direction of 
y-polarizcd spins. 



jy( x ) = s"£v x py(x). 



(5) 



where S = bc/2 is the semiconductor film cross-section 
and v x is the electron velocity in x-direction |l6j . Finally, 
Eq. J3J can be transformed to 



£ so = 7 J o dxjy(x)—. 



Here the coupling constant 7 is given by 



7 = 7o / b X iv,z = 0) dy , 



(6) 



(7) 
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where 70 = 2m*a/hS. 

From the last equation, it is seen that the spin- 
polarized flow imposes a distributed torque on the rod. In 
order to study this effect in detail we will neglect, for sim- 
plicity, the difference between elastic constants of semi- 
conductor and metal parts of NMS. As such, the equation 
of motion for the torsion angle can be then written as 



dt 2 dx 2 dx 



0. 



(8) 



where rj(x) denotes the Heaviside function, K stands for 
the torsion rigidity, and I is the moment of inertia. It is 
easy to figure out that the torque imposed by the SOI on 
NMS can be expressed as 



T = 1 / dxJ y (x)=jJ y , 
L Jo 



(9) 



and, for the S-M rod clamped on both ends, the torsion 
angle at x = L 



9l 



L(L t — L)T 
L t K 



(10) 



where L t is the total length of the rod. From Eq. © 
one can easily see that if the semiconductor film covers 
the entire length of the beam [L = L t ) and the spin 
current is homogeneous along it, the last term in Eq. (JSJ) 
turns to 0. Consequently, for a doubly clamped beam 
the solution of Eq. © is 8(x) = 0. In this case, in order 
to obtain the finite torsion angle, the NMS must include 
films with different spin-orbit coupling parameters 7, as 
in Fig. 1 where 7 = in the metal film. On the other 
hand, if J y depends on x, as in the example considered 
below, the metal film is not so necessary. In this example 
it is shown, however, that such a film can be useful as 
a reservoir for fast spin relaxation, enhancing thus the 
diffusive spin current flow through the beam. 

In order to evaluate the torque, let us adopt the fol- 
lowing simple model, which is also convenient for an ex- 
perimental realization. Namely, we assume that the spin 
current is due to spin diffusion from the left contact . The 
spin polarization P y (0) can be created there by various 
methods ranging from absorption of circularly polarized 
light to injection from a ferromagnet |3j. One more pos- 
sibility is the electric spin orientation For the steady 
state the diffusion equation reads 



D d 2 P y P y _ 
% dx 2 Ti 



(11) 



where Di and Tj are diffusion coefficients and spin relax- 
ation times, with the subscript i indicating the physical 
quantities in semiconductor (0 < x < L) (i=S) or metal 
(x > L) (i=M) regions. At the semiconductor-metal in- 
terface the diffusion current and magnetization P y /-/Vj(O) 
must be continuous, where ATj(0) is the semiconductor or 



metal density of states at the Fermi energy |17|. We 
will assume that the length of the metal part of the rod 
is larger than the spin diffusion length Im = \JP>m t ^i- 
Therefore, the spin current passes through the semicon- 
ductor film and further decays within the metal film. Ob- 
viously, in the considered example there is no charge cur- 
rent through the system. Solving the diffusion equation 
for Is 3> L and {(tmL)/(ctsIm) 3> 1, where cm and as 
are the 3D conductivities of metal and semiconductor, 
respectively, we obtain 



J y 



D s P y {0)S 



(12) 



Since the ratio <jm/<ts is very big, Eq. I|12|) is valid in a 
broad range of not very small L. 

For a numerical evaluation of the spin-orbit torsion ef- 
fect we take b = 400 nm and c = 200 nm. The SOI 
coupling constant a/h = 4 x 10 5 m/sec in GaAs [lij . 
From Eq. and Ref. [3, it is easily to obtain the spin- 
current-torsion coupling parameter 7 = 70^2^, where 
ki is a numerical factor depending on the ratio c/b. At 
b/c = 2 the factor &2 = 0.03. For such numerical param- 
eters we find 7 = 2.4 x 10~ 32 J sec. It is interesting to 
compare the torsion effect from the strain induced SOI 
with that produced by spin flips at the FM-NM inter- 
face 0). In the latter case T = h,I s , where I s is of the 
order of the spin current injected at the FM-NM contact 
when the electric current passes through it. Compar- 
ing this expression with Eq. |J3J, it is seen that at the 
same spin currents the SOI effect is much stronger, by 
the factor j/h ~ 2.2 x 10 2 . On the other hand, in H3 
the FM-NM contact can be fabricated from all metallic 
components, while our device must contain the narrow 
gap semiconductor film. In the former case NMS is able 
to carry much larger spin current, due to the weaker, by 
the factor ~ <7s/<tm, Joule heating effect. However, the 
measurement setup suggested by Mohanty et. al [12| al- 
lows to measure torsion effects produced by quite weak 
currents. For example, at eJ y = 10 _8 Amp the torque 
T = 1.5 x 10~ 21 N m, which is within the sensitivity 
claimed in Moreover, the measurement sensitivity 

can be enhanced 0] . Within our model we can evaluate 
the spin polarization P y (0) which can produce a measur- 
able effect on NMS. From Eq. JEJ, taking L = 2 /zm, the 
typical low temperature diffusion constant 300 cm 2 /sec 
and n = 10 17 cm" 3 , one obtains eJ y = 10 i (P y (0)/n) nA. 
Hence, a measurable 10 nA spin current in NMS can be 
created by diffusion of spin polarization from an adjacent 
reservoir containing only_0.1% of spin-polarized carriers. 
Various methods 0, S 13 are a ble to provide such and 
even much larger spin polarization. Higher spin currents 
are, however, restricted by the heating effects, which de- 
pend on the practical design of NMS. 

It should be noted that the torsion measurement 
method of Ref. 01 applies to a time-dependent torque 
in resonance with a NMS oscillation. For such a mea- 
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surcmcnt the spin current could be modulated in time 
by a narrow gate between the left contact and the rod, 
or by varying the spin polarization in the left reservoir, 
for example, if it is created by absorption of circularly 
polarized light with modulated intensity. 

The static torsion angle at x — L can be found from 
Eq. (|10fl . On the other hand, the maximum torsion ef- 
fect is obtained for the time-dependent spin current in 
resonance with the NMS fundamental oscillation. In this 
case, the torsion angle Or, in Eq. IjlOjl must be multiplied 
by Q/2, where Q is the resonance quality factor, which 
can be quite large in NMS. To observe this torsion angle 
it must be much larger than the mean amplitude of its 

thermal fluctuations W 86\. For a doubly clamped rod 




For a rectangular cross-section with b/c — 2, the torsion 
rigidity K = 0M7fib 3 c 0, where ft, ~ 3.3 x 10 10 N/m 2 
in GaAs material. Taking L t = 5 fim and all other pa- 
rameters the same as in the previous paragraph, Q = 10 4 
and T = 100 mK wc obtain the ratio 56 L /6 L ~4x 10~ 2 
at e.J y = 10 nA. 

We have considered a simple example of the spin-orbit 
torque effect produced by spin flux in a diffusive 3D semi- 
conductor film. It would be interesting to study other 
systems, for example, a superlattice of remotely doped 
high mobility quantum wells in the ballistic regime (L is 
less than the elastic mean free path). In such a system en- 
ergy dissipation within the semiconductor film is reduced 
and, apparently, larger spin currents arc allowable. 

In summary, we propose a nanomechanical system 
where due to the strain-induced spin-orbit interaction 
the electron spin degree of freedom can couple to NMS 
mechanical motions. We have shown that this coupling 
is strong enough to induce the measurable torsion in 
NMS when the spin polarization flows through the 
suspended nanobeam. Besides a potential for other 
possible applications, such NMS can be employed as a 
sensitive detector of spin currents and spin polarizations. 
The basic structure can be further modified to create 
devices for eventual use in spintronics as well as spin 
information processing. 
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